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FOREWORD 

This technical report documents the results of work performed by Northrop 
Services, Inc., Huntsville, Alabama, under Contract NAS8-33197 for the National 
Aeronautics and Space Administration, George C. Marshall Space Flight Center. 
The Contracting Officer's Representative (COR) for this study was Mr. A. W. 
Deaton of the Systems Analysis and Integration Laboratory, EL23. 
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ABSTRACT 

The Inertial Upper Stage (lUS) will be guided by the Boeing Company’s 
Gamma Guidance Algorithm. This report presents the results of an independent 
assessment of the algorithm’s performance in satisfying the NASA missions’ 
targeting objectives. The report includes an introduction to the algorithm, 
the results of a launch window analysis for a Galileo mission, and suggested 
improvements . 
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Section I 


INTRODUCTION 


1.1 GENERAL 

The Boeing Company is developing an onboard adaptive guidance algorithm 
for the Inertial Upper Stage (lUS). This algorithm is called "Gamma Guidance", 
and it is designed to provide real-time guidance that will satisfy all DOD and 
NASA mission targeting objectives. Ganna Guide employs a sophisticated iterative 
approach to the guidance problem^and it uses numerical onboard integration for 
some nonlinear differential equations and closed form approximations for other 
nonlinear differential equations. It supplies ignition times and thrust 
attitude for the main stages' solid rocket motors (SRM), and ignition time, 
thrust attitude and termination times for corrective translational burns of 
the reaction control system (RCS) . 

The MSA planetary missions represent a different class of lUS missions 
when compared to DOD Earth Orbital Missions. Planetary missions require earth 
escape which results in hyperbolic orbits; thus, the targeting logic and 
convergence properties of the algorithm will be different than for the earth 
orbital elliptical trajectories. Gamma Guide has provisions for solving these 
classes of problems. The Gamma Guide algorithm has the following salient 
features. 

1. It solves the two point boundary problem described by 

0 An initial position (moving with vehicle) 

o A set of selectable constraints that include targeting end 
conditions 

o A set of selectable controls encompassing ignition, attitude, 
and thrust profiles 

o A set of nonlinear equations to propagate the equations of motion. 

2. It generates ignitions times and steering angles for the Sol''d Rocket 
Motors (SRM) ; 

3. It generates ignition times, steering angles and termination times 
for the Reaction Control System (RCS) ; 

4. It is a real time, iterative, baseline guidance system that can 
satisfactorily solve the unique problem of using fixed Impulse SRM's 
to power missions from near-earth to interplanetary orbits. 
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The elgorltto has been under development for several years. References 
1, 2, and 3 contain a comprehensive explanation of the algorithm. The 
referenced explanation includes derivation, approximations, and the ass\;unptions 
used by or in the algorithm. Reference 4 contains selected portions of 
Reference 3 to familiarize the reader with the basic concepts of the guidance 
algorithm. Reference 4 also contains a listing of the double precision 
developmental code as modified for MSFC and the results of previous tests of 
the algorithm for several interplanetary missions and a geosynchronous mission. 


The algorithm has been designed to guide many different types of missions. 
It is very flexible and the Boeing B-S code for the onboard computer reflects 
this flexibility. The B-5 code, however, requires a driver code to simulate 
the inputs from the vehicle and the real world enviroiment. The code is being 
revised as necessary and the current revision Is revision I. A source program 
for this revision was obtained from Boeing, and was adapted to the HSFC UNIVAC 
1100 computer systems by MSI. 

A driver program for the B-5 code is not available, however as part of 
the B-5 simulation, Boeing supplied a six degree of freedom stand alone 
simulation (STARS) for attitude control of the lUS, that uses information 
supplied from the B-5 simulation. This code has been adapted by MSI to the 
MSFC UMIVAC 1100 system also. 

The B-5 and the STARS code were adapted to the UNIVAC 1100 system, 
although this task was not an original part of the contract. This was done 
to replace Task VIII of the contract which required data not yet available 
for delivery by the government. The change was agreed to by the COR. 

1.2 OBJECTIVES AND ACCOMPLISHMENTS 

There were nine tasks assigned to this study. These tasks were a 
continuation of the basic tasks previously investigated in Reference 4, and 
are listed below with a brief discussion of their status. 
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1.2. 1 Task 1. UPdatt Ttttbtd ProiraMt 
Raquir^Miic: 

Haka allowancea for all discrete and analog Inputs/outputs to/frcm the 
guidance algorittm that are required to aasess the perfomance of the algoritha 
from a basic guidance standpoint in the guidance sinulation progrwB. Include 
the capability of siaulating off-noninal variations in nainstage and RCS 
propulsion paraneters, lUS mass properties, end initial position and velocity. 
Also, include the capability to sioiulate mainstage ignition failures. Variations 
in navigation paruieters need not be simulated in this study; the measured 
position and velocity of the lUS can be set equal to the actual position ax»l 
velocity. 

Status: 

The testbed program has been continually updated, and Includes all of 
the detailed items required. During the course of this study, the code was 
changed to single precision calculations to simulate the on board computer. 

1.2.2 Task 2. Maintain Testbed Guidance Aieorithm 

Requirement: 

Maintain an up to date guidance algorithm consistent with the Boeing 
Aerospace Corporation latest changes. 

Status: 

NSl made one trip to Seattle to discuss with the Boeing representatives 
their latest changes to the algorithm. Additionally, NSI has kept in telephone 
coraminicatlon and discussed with Boeing various suggestions for improvements 
to the code. The MSEC IDS team has been kept Informed of the latest version. 

The improvements are presented in Section IV of this report. 

1.2.3 Task 3. Identify Mission 

Requirement : 

Identify a group of off-nominal situations for each NASA lUS mission. 

For each situation, specify the variations in the propulsion, mass, and initial 
state paraiBeters, and specify mainstage Ignition failures. Each group of off- 
nominal situations should be representative of the anticipated three-sigma 
error volume. Form a covariance matri;: using these error sources for each 
design reference mission. 
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Th« HASA mlMloiM have been in a atata of changa. The study was applied 
to the revised Galileo, ai^ to the Venus Orbiter missions, ^hasis was placed 
on the Galileo mission as a representative of the basic NASA mission objectives. 
As part of this task, NSI investigated the preload data needed for the three 
on orbit launch opportunities for the first launch window of the Galileo 
mission. The results of this study were discussed in Section 3. The covariance 
matrix source program has been updated as necessary. 


1.2.4 Task 4. Generate Guasioptimum Reference Trajectories 

KequirMMnt: 

Use SAtffiO to generate quasioptimtaa reference trajectories for the nominal 
and off-nominal situations identified in Task 3, for each mission in the 
representative set of NASA lUS missions. These reference trajectories will 
provide a basis for the analyses of the gamma guided trajectories generated 
in Task 5. 


Status: 

SAMBO was used to generate the reference trajectories for the missions 
in Task 3. Changes were made to the SAMBO code to calculate and output additional 
data needed for the MSFC Gama Guide Code. The SAMBO changes arc discussed 
in Section VI of this report. 


1.2.5 Task 5. Analyze 

Requlrenent! 

Analyze Che gaisBa guided trajectories from Task 3 and record the 
targeting accuracies, RCS propellant usages, and any abnormalities in trajectory 
profiles or attitude histories as compared to the SAMBO trajectories generated 
in Task A. Identify those trajectories that do not satisfy the targeting 
accuracy requirements, RCS propellant budgets, basic sequencing and timing 
requirements, or attitude-rate constraints. 


Status: 

This task was modified because the NASA missions do not use RCS burns 
for corrective action. The trajectories were in general well within the 
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targeting accuracla*. Soam problMui %r«r« found in basic saquancing. Thasa 
aara discuasad with Boaing and have baan rectified by the addition of a time- 
line checking saquanea in the MSPC coda. 

1.2.6 Task 6. Patch and Reaair as Necessary 

Raqulrenant : 

For each of the unsatisfactory gansBa guided trajectories (if any) identi- 
fied in Task 5, define and incorporate modifications to the gamma guidance 
logic (Fortran version) that will produce a satisfactory trajectory in the 
sami situation. Verify that the modifications are sufficient by using the 
sinmlator testbed (with the modified guidance logic) to generate a satisfactory 
trajectory. 


Status: 

All problems encountered were eliminated by adding appropriate coding. 
These changes were discussed with Boeing for possible inclusion in their 
B-S Fortran version of the on board coc- . These changes are presented also 
in Section IV. 


1.2.7 Task 7. Document 

Requirement: 

Document the results of the study in a final report. In addition to 
a general assessment of the Gatmta Guidance algorithm, doewnent the results 
on a miss ion-by-mis Sion basis. Include all suggested modifications to the 
guidance algorithm in a separate section of the report. 

Status: 

This report is in response to this task. In addition Reference 5 was 
used to document the results of a study of the preload date of the PL13A Galileo 
probe . 

1.2.8 Task 8. Modify the 6-D Protram (lUS-SIXD) 

Requlrwsent: 

Modify the six degrees of freedom spin stage computer simulation program 
built during the first contract to Include the RCS spin-up control system. 
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Ch«^ char«ce«rlstirB of the hydrazine blowdown syatm, the gyro data sampling 
and proceaaing charaeterlatica, and the open loop part of the apin~up after 
syatM aaturatlon. Documentation of the software and program users manual 
shall be In compliance with FIPS Publication 38. 

Status: 

The six degrees of freedom spin stage simulation program SIXD of 
Reference 10 was nwdified to include the basis for an RCS spin up system. 
However, the data required for completing this task was and is not available. 

To replace this task, the COR requested and NSI added the adaptations of the 
Boeing B«5 Revision I Fortran version of the on board code to the MSr'C UNIVAC 
1100 computer systma. 

The B-5 code requires a driver code to simulate the inputs from the vehicle 
and the real world environment. The driver program is not availabJe, however 
as part of the B-5 simulation, Boeing supplied a six degree of freedom stand 
alone simulation (STARS) for attitude control of the TUS, This simulation 
was also adapted by NSI to the MSFC UNIVAC 1100 system. Although the B-5 driver 
is not available, the progracmiing can be verified by Judicious use of the 
HSFC developmental code. 

1.2.9 Task 9. Provide Mission Analysis Support 

Requirraent : 

Provide a limited (no more than 20% uf the manpower) amount of analysis 
to support the MSFC/JSC lUS mission planning activities as directed by the 
MSFC contracting officer's representative. 

Status: 

NSI accompanied the COR to JSC Houston. NSI delivered and discussed the 
latest version of the MSFC Gamma Guide Development code with JSC lUS personnel. 

1.3 DEVELOPMENT CODE PHASESt STAGES AND TIMELINE 

In continuing the assessment of the algorithm, the developmental code 
of Reference 4 was used. This double pre Irion code has been obtained from 
Boeing in 19V8 and has been modified and improved by NSI. This development 


1-6 


m-L9S5 




Eng Inttrlng 4 Ttctinofogy 


Ccnttr 


code simulates the vehicle and real world input to the algorithm and is used 
as a teat bed for changes to the algorithm. During this study, the code was 
changed to single precision, so that a closer simulation to the ultimate onboard 
version would be obtained. 


The MSFC code has provision for six phases plus a terminal coast. Each 
phase consists of a coast arc followed by a burn arc for a total of 13 stages 
when the terminal coast is included. These are shown in Figure 1-1. Also 
included in Figure 1-1 is the relation of the event timeline to the Gamma 
Guide controls ^ which are shown in Figure 1-2. The development code is not 
as flexible as the B-5 code in that the timeline must reflect an RCS burn 
between two SRM burns. This is not fatal, however, since a zero second burn 
is allowed. 

A simplified flow chart of the MSFC cods is given in Figure 1-3. This 
figure also Incorporates several improv^ents suggested in Section IV o' this 
report. 

1.4 REPORT SECTION DESCRIPTIONS 

Section II of this report is a brief discussion of the type of data net ed 
by the MSFC Gamma Guide simulation and the differences between the data requ. ed 
by the MSFC and Boeing B-5 codes. 

Section III is a report of the data used in the launch window analysis 
of the Galileo mission PL13A. A preliminary report of the analysis was con- 
tained in Reference 5. Data from References 6 and 7 were also used in the 
analysis . 

Section IV is composed of suggested improvements to the algorithm. These 
have been discussed with the Boeing Company's representative involved in 
coding their B-5 simulation. 

Section V discusses the availability of source listings and cards for the 
various codes used in this study. 


1-7 


TR-1995 



£flgifl«tr<nf 4 Ttehnology Center 


Section VI is a report of changes to the computer code SAMBO of Reference 
8. These changes were made to facilitate the generating of reference 
trajectories. 


Section Vll contains the references list. 
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Figure 1-3. GAMMA GUIDE SIMPLIFIED FLOW 
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Section II 


PRELOAD DATA 


2.1 DEVELOPMENT VERSUS ONBOARD CODE 

There are tvro versions of the GAMMA GUIDE algorithm code. These ar ei 
the developmental code modified for use on the Unlvac 1108 computer system 
at MSEC, and the B-5 Revision I code being developed by Boeing for the onboard 
computer. In general, the preload data is the same for both codes, however, 
certain differences in data form and usage exists. These differences are 
shown in Table I. The remainder of this section applif’s to the data input 
to the MSEC code. 


Table 1 


BASIC DIFFERENCES INPUT DATA 


DATA 

MSFC CODE 

B-5 CODE 

At for impulse 
placement 

Integrated from 
thrust table input 

Calculated from 
input polynomial 

Target Conic 

Inertial coordinates or 
conic parameter input 
without retargeting 

Retargeting using 

polynominal 

input 

Real World, 

Gravity, Thrust, etc. 

Input for simulation 

Sensed through 
NAV system 

Times 

Referenced to program 
time zero 

Referenced to GMT 
and/or nominal 
launch time 

Control Lambdas 

Discrete data points, 
referenced to inertial 
system at time zero. 

No retargeting 

Retargeting with 
polynominal s 
referenced to 
nominal launch 
times and on-orbit 
delays. 


2.2 MSFC SIMULATION 

A complete description of the input data to the MSFC code can be found 
in Reference 4. Additionally Reference 5 contains a discussion of targeting 
preset data for the Galileo mission. Portions of References 4 and 5 are 
Included here as an introduction. A broad overview of the preload data required 
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by the MSFC GAMMA GUIDE simulation program is shovm in Table 2. The data is 
divided into 3 categories. Types I and II are the preload needed bv the 
algorithm, while Type III simulates the real world aspects for supplying 
navigational data for the algorithm. 

Table 2 

INPUT DATA TYPES 

I. NOMINAL NON-CONTROL PRELOAD 

Initial position 
Target orbital data 
Various control flags 
Predicted SRM thrust tables 

" SRM burn time 

" RCS ISP 

" RCS thrust equation constants 

" Drop weights 

" Stage weights 

Used "as is" by the algoritiwi with (type II below) to predict 
future states. 

II. NOMINAL CONTROLS (Lambda's) 

Impulsive delta velocities 
Impulse times 

Modified by the algorithm and used with (type I above) to predict 
future states. 

III. NAVIGATION DATA FOR SIMULATION 

SRM thrust tables 
SRM ISP's 
SRM burn times 
RCS ISP 

RCS thrust equation constants 
Used by the navigation system with the computed 
controls of (type II above) to update the actual state. 

A listing of a sample case input for the MSFC code is contained in 
Appendix A. 
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Section III 

PL13A GALILEO PROBE 


3.1 MISSION PROFILE 

The following mission scenario and ground rules are extracted from 
Reference 6. 

Begin Quotation: 

A. Mission Scenario 

Before going into the details of the launch window trajectory generation, 
it would be helpful to summarize the salient features of the Space Transporta- 
tion System (STS) and the Galileo mission. The lUS/Galileo spacecraft will 
be carried to a deployment orbit by the Space Shuttle. The first two burns of 
the Space Shuttle Orbiter On-orbit Maneuvering System (OMS) will place the 
Orbiter /cargo Into an 100 x 100 nml orbit. After 3 to 4 revolutions in the 
100 X 100 nmi orbit, a third OMS burn will place the Orblter/cargo into a 
100 X 250 nml deployment orbit. After deployment of the lUS/SC, the lUS, 
which is made up of two large solid propellant motors (21,000 lbs. class) 
and one small spin stabilized solid propellant motor (6,000 lbs. class), will 
execute the 1st large motor burn across perigee of the deployment orbit and 
thereby place the lUS/SC into a 100 x 3200 nmi elliptical transfer orbit. After 
approximately one revolution the 2nd large motor will ignite and burn to cut- 
off. The 3rd solid propellant motor will ignite upon completion of a 270 
second coast period. The 270 second coast period will be varied to compensate 
for over /under performing motors. The 3rd stage is spin stabilized to 70 RPM 
with the Reaction Control System (RCS) . The 2nd stage and RCS/avionics section 
is staged prior to 3rd stage ignition for AV efficiency. Upon completing the 
3rd stage burn, the spacecraft will despin the configuration. Approximately 
10 days later the S/C will execute a AV course correction to enable the S/C 
to swing-by the planet Mars on the way to the planet Jupiter. Upon arriving 
at Jupiter (approximately 31/2 years later), the Galileo spacecraft will 
brake into orbit about Jupiter and encounter several of the Jovian satellites. 

B. Space Shuttle Orbiter/IUS/SC Groundrules and Assumptions 

The integrated Orbiter/IUS trajectories for the January 4, 1982, to 
January 14, 1982, Galileo mission launch period were generated under the following 
groundrules and assumptions for the Orbiter, lUS, and S/C; 
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ORBITER: 

1. Approximate Midnight Launch (1/4/82 - 1/14/82) 

2. The daily launch window should be defined as that period of time 
determined by the launch azimuth constraints of 85*£ AZ <,95* . 

3. No yaw steering requirements shall be Imposed by the planetary 
targeting objectives. However, any real time launch environment or 
vehicle perturbation that causes a drift out of a plane will have to 
be corrected by Space Shuttle yaw steering. 

4. The Space Shuttle flight profile shall use a two (MS burns ascent 

to a 100 X 100 nml orbit followed by a minimum coast time of 4.5 hours 
prior to the third (MS burn that achieves the 100 x 250 nmi lUS/SC 
deployment orbit. 

5. The minimum coast between Orbiter insertion into the 100 x 250 rani 
deployment orbit I’.nd lUS/SC deployment shall be 3 REVS in orbit. 

6. The Orbiter shall provide the lUS initialization data (time, position, 
velocity, and IMU alignment angles) to the lUS computer prior to 
deployment . 

7. The Orbiter shall provide pointing during the coast phase prior to 
deployment consistent with the thermal conditioning requirements 
of the lUS/SC. 

8. The Orbiter shall provide the translational AV for the lUS/SC separa~ 
tlon after deployment and prior to the lUS first stage Ignition. 


lUS: 

1. The lUS first stage propellant off-load shall be determined from the 
following conditions: 

a. Maximum C3 of the launch period. 

b. The third on-orbit launch opportunity. 

c. 30 second on-orbit ignition delay of the 
first SRM. 

2. Reference trajectories shall be developed for the 1st day of the launch 
period, fifth day, and last day. 

3. The tracking analysis shall be based on the reference trajectories 
from the fifth day, launch azimuths of 85®, 90®, and 95* and the 
first, second, and third on-orbit ^aunch opportunities. 

4. The guidance (steering) analysis shrll be based on the minimum C3 
launch period (energy management exercised) , and the third on-orbit 
launch opportunity. 

5. The navigation analysis shall be based on the minimum C3 launch 
period, and the 3rd on-orbit launch opportunity. 

6. The spin stage accuracy analysis shall be based on the minimum C3 
launch period, and the 3rd on-orbit launch opportunity. 
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7. Th« lUS/SC daployaant state shall reflect the Space Shuttle ascent pro- 
file to a 100 X 100 i»l orbit and subsequent 100 x 250 nni deployment 
orbit. 


GALILEO SPACECRAFT: 

1. 2060 kilogram spacecraft 

2. Project Galileo Targeting Specifications (625-120 dated 9/15/78). 

Ei^ Quotation. 

This mission starts from a low parking orbit, with an SRM burn at perigee 
to transfer to a higher energy elliptical orbit, coasts backs to perigee, then 
has a second SRM burn followed by a coast, and then the third SRM burn which 
injects into a hyperbolic orbit that goes to Jupiter. The mission is shown 
in Figure 3-1 and further described in Reference 6. 

The launch azimuths for this mission are restricted to the range from 
85 to 95 degrees and there are three on-orbit launch opportunities for each 
window. The first opportunity is a dog leg to the left maneuver, occurring 
approximately one orbit before the nominal in plane opportunity for the daily 
launch window. The third opportunity is a dog leg to the right maneuver 
occurring approximately one orbit after the nominal. The launch window and 
opportunity selected will determine the preload values needed for the slim.tlation. 
Most of the data described in Section II remains fixed. The variable data 
include: 

o The initial state (position, velocity), 

o The predicted stage weights, 

o The target orbital data (desired mission constraint values) , 

o The time at which the mission constraints are evaluated, 

o The nominal controls (Lambdas) which Include the impulsive delta velocity 
and Impulse time for each SRM burn (RCS also if used). 

The B-S onboard code will have provisions for updating these as a function 
of launch time and on-orbit delay. The study described in Reference 5 was 
undertaken to test the feasibility of the retargeting procedures described 
in Referenc-' 7. 
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3.2 6ENERATI0N OF PRELOAD DATA 

The prograa SAIOO fsimplex Algorithm for Multiple Burp Optlnizatlon, 
Reference 8) wee need, with date frra Reference 6, to genernce reference 
trejectories for the Galileo Hiseion PL13A. Using azimuths of 85.0, 87.5, 

90.0, 92.5, and 95.0 degrees, references were obtained for the three 
opportunities of the first launch window. The date from these runs •^ore then 
used sa trial data load for GAMUA GUIDE. These 13 cates all ran satisfactorily 
using at aost four interations to convarge on the first <pre SEMI) solution. 

The converged data is shown in Tables 3, 4, 5, and 6. This data can be divided 
into: 

0 Predicted Stage Weights 
o Target Conic Elements 
0 Control Values 
0 Initial State 


3.2.1 Predicted Staec Weights 

These values are shown in Table 3 and are used to predict AV output from 
each SRM. The AV from SRMl may be input as DELVl. If not input it will be 
calculated Internally by 

AVI - Go Isp^ In (^) 

The values of AV2 and AV3 are obtained, in the MSFC version, from the predicted 
thrust Cables for the SRM. In this study, it was found Chat the values of the 
predicted stage weights and AVI for opportunity 2 could be used as estimates 
for opportunities 1 and 3 also. Note that this implies that these values 
can be fixed or set as a function of the window only. 


Table 3. PREDICTED STAGE WEIGHTS AND AVI 



1 

OPPORTUNITY 

2 

3 

AVI 

3390.09 

3391 .00 

3391 .89 

HT3 

55497.41 

55485.72 

55473.60 

WTX 

38774.55 

38762.56 

38750.43 

WT5 

36502.46 

36490.47 

36478.35 

WT9 

11347.48 

11347.48 

11347.48 
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3.2.2 Tarftt Conic Elements 

The target conic el^ente are referenced to the inertial coordinate 
•ystra ^ich is eatabllshed at program time sero. The elesMnts are defined 
in Figure 3-2 and their values are given in Tables 4. 5 and 6. They are: 


0 

a 

Sraimajor axis , 

0 

e 

Eccentricity , 

0 

1 

Inclination , 

0 

w 

Argxaaent of Perigee . 

0 

n 

Ascending Nodal ^gle 

o 

n 

True Anomaly, 

0 

TIME 

Time . 


It should bt; noted that these variables are almost all constant as a function 
of the azimuth. The exception is the ascending nodal angle which varies as a 
function of time and the earth's rotation rate. This is a result of the 
inertial coordinate system being fixed at time zero. The target orbital data 
was therefore revised to the averages shown in Table 7. 
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Rj IDS Position Vector (Navigation System) 

X.Y.Z Basic Earth Centered Inertial (ECI) 

Coordinate System For Guidance Computations 
i Inclindtion of the Target Orbit 

P Perigee of the Target Orbit 
w Argiment of Perigee 

n Right Ascension of the Ascending Node 

n True Ananaly--Vehicle Position with Respect to Perigee 


Figure 3-2. THE GUIDANCE COORDINATE SYSTEM AND THE TARGET 
ORBIT DEFINITION 
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Table 7 


REVISED TARGET ORBITAL DATA 

Sanimajor axis a -20123482 
Eccentricity e 2.08879 
Inclination 1 28.59204 
Argument of Perigee u 71.07624 
True Anomaly n 93.88088 


ASCENDING NODE ANGLE 


AZIMUTH 

OPPl 

0PP2 

0PP3 

85.0 

-83.91125 

-83.98913 

-83.87525 

87.5 

-89.22364 

-89.20139 

-89.17751 

90.0 

-94.15586 

-94.13065 

-94.10520 

92.5 

-99.72639 

-99.70639 

-99.67873 

95.0 

-104.47325 

-104.45130 

-104.37535 


3-11 


Eiif Inter ifif 4 Ttehnology Center 


TR-1995 


3.2.3 Control Lambdas 

The converged values of the control variables are shown In Tables 4, 5, 
and 6. The RCS values of zero are needed by the MSFC code, but are not used 
by the B-5 code. 

The data for the SRMl burn were then modified to reflect a uniform 
AVI of 3391 ft/sec corresponding to the value for Opportunity 2. Uniform 
times for the RCS burns were also set. The reduced data are shown in 
Tables 8, 9, and 10. 

3.2.4 Initial State 

The initial state values for the MSFC simulation are as shown in Table 
11. These values were used as functions of the launch azimuths. 
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Table 11. 

INITIAL STATE 






AZIMUTH 




85.0 

87.5 

90.0 

92.5 

95.0 

LAT 

28.63360 

28.249024 

28.40993 

28.33325 

28.495628 

LON 

14.84890 

11.85565 

1.87490 

1.28092 

358.76512 

ALT 

617800.6 

622245.4 

618443.9 

622255.58 

626278.21 

AZM 

92.95830 

94.02688 

91.62042 

93.98590 

95.12122 

GAMMA 

.43980 

.46922 

.44485 

.46896 

.49189 

VEL 

25815.41 

.2581035 

25814.70 

25810.31 

25805.58 


3.3 INTERPOLATION FOR OTHER LAUNCH AZIMUTHS 

Linear interpolation was used on the reduced data to supply initial data 
load values for the launch azinuths of 86.25, 88.75, 91.25 and 93.75 degrees 
for each of the three opportunities fk«r the first window. These cases ran 
satisfactorily with initial convergence obtained in not more than six iterations. 

3.4 RETARGETING 

The proposed method for resetting the GAMMA GUIDE preload data as 
necessary to compensate for launch, deployment and on-orbit delays is given 
in Section 7.1 of Reference 7. The on-board retargeting consists of four 
primary tasks: 

o Updating the SRM Impulse directions, 
o Updating the SRM Impulse time, 
o Updating the mission constraints, 
o Updating the initial state. 

Only the first two of these tasks are addressed in Reference 7. The last 
task is assumed to be taken care of by the onboard navigation system. For 
the purpose of retargeting, the SRM impulse directions and times are calculated 
in a space centered inertial system (such as the M-50) and then converted over 
to thfe ECl system. This takes advantage of the fact that the Impulse directions 
can be specified as a yaw angle a followed by a pitch angle B as shown in 
Figure 3-3. Thus if the spin axis of the two coordinate systems are aligned 
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AV * (AX.iV.AZ) 

Aii « aV Cos 6 Cos a 
aY = aV Cos g Sin a 
aZ « aV Sin g 


Figure 3-3. IMPULSE DIRECTION FROM a.g 


ingiitetfittg 4 TvcMsIogy Ctnttr 


TR-1995 


the transfornietlonH: 
. * 


a 

8 




can be used to relate a end 8 in the ECl systen to o', 6* in the retargeting 
coordinate system; «rhere t is a suitable time relating the two systras. In a 
siAilar manner, impulse times can be reset 

TS»f • TSSM' + t. 


From this preliminary study, it appears that the preload data can be 
impl«Bcnted either in tables with linear Interpolation or we can use linear 
equations. Although the study was accomplished in terms of launch azimuth 
and opportunities we note that the driving independent variable is really the 
time at which the Earth Centered Coordinate Syst«8, ECI, Is initialized. This 
is particularly true for the ascending node angle Q and the firing angle a 
for each burn. Note that: 

o The target elements a, e, 1, u and n can be fixed. 

o The target element Q (ascending node angle) should be a function of 

the time at which ECI is initialized. 

o The Initial values of a and B for each burn can be set to the values 

for opportunicy 2. These are close enough for the algorithm to correct 

for opportunities 1 and 3. 

0 The values of AVI, PWT3, PWT5, W9 and WTX can be fixed at the values 
for opportunity 2 also. 

o The on->board retargeting for lUS GAMMA GUIDANCE as given in Section 
7.1 of Reference 7 is feasible for the Galileo Mission PL13A. 


3.5 PRECISION OF CALCULATIONS 

This data preload investigation used the double precision version of the 
MSFC GAMMA GUIDE code. Sample cases using a single precision code have been 
also run satisfactorily. In single precision, minor change in impulse times 
and directions were encountered and the cases ran in less computer time. The 
MSFC Univac 1108 uses a 36 bit computer word while the onboard computer is 
expected to use only 32 bits. This is only expected to affect those operations 
Involving repeated divisions or the square root. 
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3.6 RELIABILITY AND NASA NISSIONS' OBJECTIVES 

Rafsrmec 4 pr«s«nt^ the result* of s study by NSI investigating t'ls 
algoritlsi. The report detailed various perturbation in the input parwaeters 
for the code, nie perturbations were at least as severe as the expected 3 
sigSM value*. This previous stwly concluded that GAIftIA GUIDE could satisfactorily 
flMSt ths NASA aissions* objectives. This conclusion has been further strengthened 
the results of this present study utilising the iaproveaents in the next 
section. 
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Section IV 

SUGGESTED IMPROVEMENTS 

4.1 INTRODUCTION 

Several i^roveaenti to the CAWIA GUIDE elgoritlm were previously 
suggested in Reference 4, These have been incorporated in the MSFC code. 

4.2 MINIMUM DELTA LAMBDA TEST 

When the program was changed to single precision coding, an instability 
sotMtiaes occurred in the calculations of the changes in the control variables 
5^ with respect to the changes in the constraints. The constraint deviations 
6Y would oscillate between a negative and a positive value just out of 
tolerance. During this time, it was noted that % was extremely small, 
indicating that the best possible solution in single precision had been 
obtained. suggested that a test for minimum be inserted into the 

algorittUB. This test is shown in Figure 1-3 and has been made by Boeing in 
the B-S version. 

4.3 MAXIMUM ITERATION RESET OP LAMBDA 

Whenever the maximum iteration limit is exceeded, generally the best 
possible solution for the control variables X has been obtained. However, 
there may arise cases in which the solution may be erroneous. To eliminate 
this probl«n, it is suggested that on beginning a guidance update cycle, the 
control vector X be saved and if the iteration limit is exceeded this stored 
value be used. This change has been made by Boeing in the B-5 version, 

4.4 TIMELINE CHECKING 

This problem has been previously discussed in Reference 4, In which a 
minlsaim coast test to Insure that an overlap In the time of end SRM2 burn and 
begin SRN3 burn did not occur. This test should be extended to the coast 
between any two burns. Coding for these tests have been added to the MSFC 
program. The timeline for the MSFC code is shown in Figure 1-1. 
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S«ctlon V 

PR06RAM AVAILABILITY 


5.1 LlSTINBSf HA6NETIC TAPE AND CARD SOURCE 

Th* prograiM /codas in Tsbla 12 wars usad in this study and have been 


stored on a single PUR tape in the MSFC Tape library. 

Se 

Appendix B 

their use. 

Table 12. PROGRAM FILES 



PROGRAM 

USE 

FILE 

CARDS 

APPROX. 

»^FC Otv. 

Oevelopment of Gamma Guide 

1 

4000 

NSI SIXD 

Spin Stage 6-0 

2 

1200 

Boeing B-5 

Onboard Gamma Guide 

3 

5000 

Boeing Stars 

Stand Alone Attitude 6-0 

4 

4000 

NSI NOB 

Thrust Tables for Gamma Guide 

5 

400 


These prograns are still under development, and are therefore changing. 

A copy of the latest listing can be obtained by contacting the Contracting 
Officer's Representative, Mr. A. W. Deaton of the System Analysis and Integra 
tion Laboratory, EL23, Marshall Space Flight Center, Alabama. 

A magnetic tape suitable for use on a Univac 1108 system or a card 
deck c^n also be supplied to authorized users. 
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Section VI 

SAMBO MODIFICATION'S 


6.1 MISSION PLANNING MODIFICATIONS 


Modifications were made to the SAMBO trajectory program of Reference 8 
to provide output data consistent with misslou planning efforts. Usual 
trajectory profile data, in terms of mean orbital elements, continue to be 
printed in the first row. However, three additional rows of output were 
added. The second row gives 


R(km) 

V(km/s) 

y(’) 

LAT(“) 
LONG(") 
ALT (km) 
SRR(km) 


radius magnitude of lUS/SC 
velocity magnitude 

flight path angle with respect to the local vertical 
earth geocentric latitude of lUS 
earth longitude of lUS 

altitude of lUS above the earth ellipsoid. 

slant range from the launch site along the 
earth surface to sub-satellite point. 


The third row gives position and velocity vectors, in terms of the M50 
Cartesian coordinate system, and the magnitude of acceleration. The M50 
coordinate system is the Aries-mean-of-1950, Cartesian coordinate system with 
the origin at the center of the earth. The epoch is the beginning of Besselian 
year 1950 or Julian ephemerls date 2433282.423357. The XM50-YM50 plane (see 
Figure 6-1) is in the mean earth's equator of epoch and the XM50 axis is 
directed toward the mean vernal equinox of epoch. The ZM50 axis is directed 
along the earth's mean rotational axis of epoch and is positive north. The YM50 
axis completes the right-handed system. The fourth row of trajectory data gives 

AZI(“) ... inertial heading angle of lUS Inertial velocity vector 
measured In the horizontal plane clockwise from north. 

g(®) ... the angle between the line from the earth's center to the 

sun and the lUS/SC's orbital plane. 


A modification was also made to the tracking analysis summary option in 
the SAMBO program. Eleven ground tracking stations and three tracking satellites 
are available. Output data printed for each acquisition period has Included 
azimuth, elevation, range and range rate of the vehicle with respect to the 
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■tationCs) selected. To this output has been added the cone angle between the 
vehicle longitudinal axis and the line-of-sight to the tracking station. 


6.2 CONVERCION OF DATA 

A user's guide to the conversion of reference data from SAMBO to input 
for Gamma Guide can be found in Appendix C. 
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Appendix A 
SAMPLE DATA LOAD 


This appendix presents a short listing of a sample data load for the 
MSFC Camilla Guide development program. Other sample data cases are stored on 
the PUR Tape described in Appendix B. 
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Appendix B 
USE OF PUR TAPE 


original page is 

OP POOR QUALITV 


The prograaw/ codes of Table 12, repeated below, have been stored on a 
single UNIVAC 1100 PUR Tape. This tape Is In the program data format of the 
1100 system and can be accessed by the assignment 


@ASG,TJ.TAPE,U9S,.XmX, . 


tdiere XXXXX Is a number that can be obtained from the Contracting Officer's 
Representative, Mr. A. W. Deaton, EL23, Marshall Space Flight Center, Huntsville, 
Alabama . 


Table 12. PROGRAM FILES 


PROGRAM 

USE 

FILE 

CARDS 

APPROX. 

MSFC Dev. 

Development of Gamma Guide 

1 

4000 

MSI SIXD 

Spin Stage 6-D 

2 

1200 

Boeing B-5 

Onboard Gamma Guide 

3 

5000 

Boeing Stars 

Stand Alone Attitude 6-0 

4 

4000 

NSI MOB 

Thrust Tables for Gamma Guide 

5 

400 


As a user aid each file has a data element called CHECKS vdtlch will execute 
a sample case. As an example, the following run stream will read In and 
execute a sample SIXD Case. 

@ASG , T J . TAPE , U9S , . XXXX . . 

SREWIND.TAPE 

OMOVE.TAPE.l 

@COPIN.TAPE,TPF$ 

0FREE,TAPE 
0ADD. CHECKS 
0FIN 
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Th« codes Bosing B-5 Boeli^ STARS both require a larger tmporary 
file, if chaise are to be made. This can be obtained by 

@FREE.TPF$ 

lASC.T.TPPS, till 1000 

Then repeating the run strean above with either 
f M0VE,TAPE,2 CFor 8-5) 

or 

I }R}VE.TAPE,3 

To obtain a listing, the MSFC utility program LISTIT can be called by 
replacing the ADD option with 

@XQT,AL.SYS$*MSFC$ .LISTIT 

Inserting the card 

^PUNCH.SYMBOLICS 

After the @XQT above will punch out the source statements. 
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A»»tndix C 

6AHMA 6UIDE INPUT DATA FRON SAMBO 

Rcf«r«nee 9 present^ • diftcustlon of converting reference date fr«B 
trejeetoriea generated uelng SAMBO* to a fora auitable for GAIMA (R}XDE. Ifoat 
of this procedure la atill valid* but HSI haa extei^ed the SAMBO capability 
and the lJ^>ulBive aolution ie now uaed by SAMBO to create a finite burn 
aolutlon. 1%ua the tiaes of beginning and e»i of each SBM burn need tyot ba 
calculate. To further ain^lify the input* the tine coluan of the Block 2 
ii^tut Bay be aet to aero ai^ the prograa will aupply theae tinea. Aa an 
additional aid* the iaq>ulBe del tea for a SlOf burn can be aet up by uaing 
option 18 in the MSFC veralon. Aa an exanple* for SKMl, the naneliat block 
ii^t 


XLO(l) - 3391.0, XL0(2) - 0.0, XL0(3) - 0.0 

with the Option 18 input below 
...18 

1-..71.2816.„..7.9603 

will aet 


XL(1) - 1077.7442 XL(2) - 3180.6942 Xi.(3) - 469.6091 


so that 


AVI - 3391.0 
a ■ 71. 2816 
8 - 7.9603 

The notes below, with thc: description of Input contained in Reference 4, 
will assist the user in converting data from the SAMBO trajectories to 
input for GAMMA GUIDE. 
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SAMIO TO GAMHA 


T1-199S 


1) Convert S— bo thmot tobloo 

o Punch M given into GAMMA Pom, Run in MOB prograa with tine-chrust 
fectore of 1. This converts to English units ss input to MOB is in 
■etrie but its output it li^ilish. 

o Select, froB output, 15-20 points for Genoa Guide and adjust thrust 

factor to obtain require total iopulss, and ISP to obtain oats overboard. 

o Iscrset TBl, TB2, TB3, ISPl, XSP2, ISP3, DELVTl, DELTl, DBLT2, DELT3 


2) Set up Initial point. Sanbo units are Mtrlc. 


Radius of Earth 
Gravity 

LAT - Latitude in degrees 
LON - Longitude in degrees 

FOR ALTITUMS (Sai^ ALT is above ellipsoidal Earth) 

RMAC/.3048 - 20925662. 

VMM;/. 3048 • Velocity ] 

FLTPA Flight Path Angle in degrees 1 

AZ AcisHith in degrees. I 


o Conversion factors 
.3048 a/ft 
.45359237 kga/lb 
20925662. ft 
32.174048 ft/sec^ 
o Initial position 


3) Set up weights 

o Convert ness at beginning and end of each burn to obtain 

PWTl 

PWT3 

PWTi 

0WSRM2 

CWSRM3 

WTX 

i 

I o These nsy be obtained fron MOB run. 


4) Set up Laobdas 

o XL (4) - Tine burn 1 end - OELTl 
XL(12) - Tine burn 2 end - DELT2 
XL(20) ■ Tine bum 3 end - DELT3 
o Set RCS psu^o burn tines 
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w 


o 8«t iapulsR daltM 

aa) - Dam xim - o.o xlo) - o.o 

tt(9) ■ mm a(io) - o.o xi(ii) - c=o 

tt<i7) - oam a(i8) * o.o xl<i9) - o.o 

o Set burns to desired elphe, beta by using option 18 in MSFC version. 


5) Set tentinel position to YB^ 

YBAK - XXXGB, YYYCE, ZZZCE, XXOCE. YYDGE, ZZZOGE 
these ere elreedy in teglish Units 

TEW ■ Finei tine 
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